A NOTE ON THE TENSOR PRODUCT OF TWO RANDOM UNITARY 

MATRICES 



TOMASZ TKOCZ 



Abstract. In this note we consider the point process of eigenvalues of the tensor product 
of two independent random unitary matrices of size m x m and n x n. When n becomes 
large, the process behaves like the superposition of m independent sine processes. When 
m and n go to infinity, we obtain the Poisson point process in the limit. 



1. Introduction 

In quantum mechanics the time evolution of two noninteracting subsystems can be 
described by an operator e**^ ® e**^ , where H and H' are Hamiltonians of the subsystems 
(see e.g. chapters 2.2 and 3.1 in [BP]). In apphcations, the unitary operator e**^, which is 
a priori comphcated, is replaced by a random unitary matrix, to make a model tractable. 
This powerful idea goes back to E. Wigner. Here by a n x n random unitary matrix we 
mean a matrix drawn according to the Haar measure on the unitary group U{n). From 
this point of view it seems natural to study asymptotic local properties of spectra of the 
tensor product ® -Bn of two independent m x m and n x n random unitary matrices, 
to which this short note is devoted. The note, in a sense, continues the investigations 
commenced in [TSKZZj . 

Some preliminaries are presented in the rest of this section, and the main result is stated. 
The proofs are provided in the next section. The last section is devoted to some concluding 
remarks concerning the tensor product of more than two matrices. 

1.1. Background and notation. For a simple point process r on M we denote its k-th 
correlation function, when it exists, by p^'' (for the definitions see e.g. |HKPVj ). Let us 
introduce three point processes 11, E, and S„. By 11 we shall denote the Poisson point 
process on M for which = 1 for all k. By S we shall denote the sine point process on 
M which has the correlation functions 

(1) p'i\xi, . . . , Xfc) = det [Q{xi, Xj)]^.^^ , 

where the sine kernel Q{x,y) = q{x — y) and q reads as follows 
. , , , sin(7ru) 

ITU 

Given a n x n random unitary matrix with eigenvalues e*^^, . . . ,e*^", where C,i € [0,2-7?) 
are eigenphases, we define the point process H,„ = {^i, . . . ,^n}- It is well known that this 
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process is determinantal with the kernel Sn{x,y) = — y), where 

. ^ 1 sin (^) 
2tt sm (^Ij 

i.e., 

(4) ...,Xk) = det [Sn{xi, Xj)]^.^^ . 

Since {^uj )■ q{u), when n becomes large, the process — vr) of the rescaled 

eigenphases of the n x n random unitary matrix locally behaves as the sine process S. 

By superposition of two simple point processes \1/ = {ipi, . . . , iPm}, $ = {01, • • • , 4>n}, 
M, N < oo we mean the union \1/ U $ = {ipi, . . . , ipM, 0i, • • • , (Pn}- 

1.2. Results. Given two independent m x m and n x n random unitary matrices A and 
A' we get two independent point processes of their eigenphases = {^i, • • • ,^m} and 
= . . . , respectively. We define the point process Em ® of the eigenphases of 
the matrix A^ A' as 

(g) = {^i + ^'j mod27r, i = 1, . . . ,m, j = 1, . . . ,n}. 

It has been recently shown |TSKZZl Theorem 1] that the process ® S^) behaves 

locally as the Poisson point process on ]R_|_. We refine this result and investigate what 
happens when n becomes large with m being fixed, or when both m and n becomes large 
but not necessarily m = n. 

Theorem. Let and be point processes of eigenphases of two independent m x m 
and n X n random unitary matrices. Let Si, ... , be independent sine processes and let 
II be a Poisson process on R. Then for each k < n the k-th correlation function of the 
process ® exists and 
/ ) (fc) (fc) 

{aj PmSiU...UmS„' 

uniformly on all compact sets in M^. 

Remark 1 (Weak convergence). According to |HKP V] . by a point process on M we mean 
a random variable with values in the metric space A^(M) of cr- finite Borel measures on M 
(counting measures correspond to locally finite subsets of M) endowed with the topology 
generated by the functions /i i— )■ / /d/x for continuous, compactly supported /. We say 
that a sequence of point processes (r„) converges in distribution to a point process r if the 
law Un of Tn couverges weakly to that of r, say u, in the space A^i(A^(M)) of probability 
measures on Jli{R), i.e. / /dz/„ — )• f fdu for any bounded continuous function on A^(M). 
Clearly, these integrals can be expressed using correlation functions, hence the theorem 
implies the convergence in distribution of the considered point processes. 

Remark 2 (Heuristic behind ([a])). In view of the mentioned theorem from |TSKZZj result 
(jb]) should not be surprising. Neither is ^ as in the simplest case m = 2 we have 

S2 ® ={^1 + ^[ mod27r, . . . , 6 + C mod27r} 

U {6 + mod27r, . . . , 6 + C mod27r}. 
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After sliifting and rescaling we end up witli two families of the rescaled eigenphases of 
a n X n random unitary matrix which differ roughly by a large shift ^(^i — ^2) which 
is independent on the matrix. That makes the families independent and in the limit, 
according to _ x > they look like sine processes. □ 

Remark 3 (Superposition of many sine processes becomes a Poisson point process). Notice 
that for any independent copies $1, . . . , of a point process $ we have 



(m — p)\ 

P=l 7re6(fc,p) ^ j = l 

where (S{k,p) is the collection of all partitions into p nonempty pairwise disjoint subsets 
of the set {1, . . . , k}. By this we mean that if vr is such a partition then tt = {tti, . . . , VTp}, 
where iiq = {n^q, 1), . . . , 7r(g, tl7i"g)} is the g-th block of the partition tt. 

Along with the fact that if we rescale, p^^(a;) becomes t^p$^ (i^)' previous obser- 
vation yields 

P=i 7re6(fc,p) ^ i=i ^ ^ 

When m goes to infinity we thus get 

lim pSiu...umE,„(a;) = lim JJp^^^ ( — {xi)ienA = 1 = p{f^ 



It retrieves the special case of a quite expected phenomenon put forward in |CD] . Namely, 
the authors say "[...] a Poisson process can be viewed as an infinite superposition of deter- 
minantal or permanental point processes" (see Theorem 4 therein and the two preceding 
paragraphs). Regarding Theorem ([a]) that implies 

lim lim „, ^^ = \. 

Note that in the second part of the theorem we establish a stronger statement, that letting 
the dimensions of two independent random unitary matrices to infinity reduces all the 
correlations in their tensor product. □ 

2. Proofs 

For the sake of convenience, let us recall a few basic facts which shall be frequently used. 
Note the following easy estimate (for the definition see ([3])) 



(6) sup 



27r 

) 

n 



1. 



Combined with Hadamard's inequality (see e.g. (3.4.6) in |AGZ] ). it allows us to bound 
the correlation functions, 

l.k/2 

(7) sup 4l{x) < k'^/^sr.lt = -—n\ 
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2.1. Proof of Theorem ([a]). Let Om,n = ^i^m ® 2^ — vr). Fix a natural number k. 
Since we will let n go to infinity, we may assume that k < n. First we show that there 
exists functions p^j^ ^ : — > [0,oo) so that for any bounded and measurable function 
/ : — ^ R we have 



where the summation is over all ordered fc-tuples {6i, . . . ,6k) of distinct points of 9m,n- 
This will prove that PqI^^ are the correlation functions of Om,n- Then we will deal with 
the limit when n ^ oo. 

Fix /. Since for each s = 1, . . . , k, 9s = ^{iis mod27r — tt) for some G {1, . . . , m}, 
js G {1, . . . , ra} we can write 

'l- mod27r — tt) 
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jeji,...,™}*^ 

where the second sum is subject to fc-tuples z, j such that the pairs (ii, ji), . . . , {ik,jk) are 
pairwise distinct. For sure it happens when all the jVs are distinct. Call these choices of 
i and j good and the rest bad. So 

i,j good i,j bad i j 

First we handle the good sum. Some i^'s may overlap and we will control it using partitions 
of the set {1, . . . ,k} into p < k Am nonempty pairwise disjoint subsets (see Remark E] for 
the notation) so that is = it whenever s and t belong to the same block of a partition. We 
have 



kAm 



E E f = E E E E E /■ 

good i,j p=l 7rGS(fe,p) distinct distinct 

V(i,i) v,V(p,i) iivj'fc 

The sums over z's and j's have been separated. Therefore taking advantage of independence 
as well as recalling definitions of the p-th and k-th correlation functions of and SJ^ we 
find 

^ I] f = J2f I /f(^(a;^(.)+Z/smod27r-7r))' \ 

• • • , a;p)pH';^(i/i, . . . , i/fc)dxi . . . dxpdyi . . . dyk, 

where we note 7r(s) = q <^=^ s G vr^. Finally, we need to address the technicality 
concerning the addition mod27r. Keeping in mind that we integrate over [0, 27r]^ and 
[0, 27^]^ we consider for 7] G {0, l}'^ the set 

[/^ = |x G [0, 27r]P, y G [0, 2Tif] \/s < k x^^s) + ys < ^tt if r]s = 0, and 

Xn(s) + ys > 2tt if r]s = 1 \. 
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Then on [/^ we have ^^^(s) + Us niod27r = x.„[s) + Vs — '^'^Vs, thus changing the variables on 
Urj SO that Zs = ^(a;7r(s) +ys- 27r?7^ - vr) we get 



good i,j \ p,7v,r] 



[0,2tt]p \ / 



where y^iz, x) = ^z^ - x^(s) + Svrr^s + vr, 

f 2-K 277 

14 = < s G M^; ys < k Zs + ^nrjs — n < Xj^ts) < Zg + 27r?7s + tt 

mn mn 

and 

Wr^ = {z e M^; Ms <k Zs < mn/2 if r/^ = 0, and Zg > -mn/2 if rjs = l] . 
Summarizing, we have just seen that the correlation function Pel„(-^) t^ikes on the form 

(8) PtiS'A = Y.^wM I V.(x)p£(x) (^)\^^{y{z^x))dx + B^,^{zl 



where the term Bm,n corresponds to the sum over bad indices ^^b^d i j f ■ -^^ the same 
kind of reasoning we show that roughly 



p=l q=l ne6{k,p) rj 
re&{k,q) 



p^){y{z, x))dx, 



where the sums are over appropriate partitions and W^, are suitable sets which appear 
after changing the variables. Now, by ([7]), 

K^) IIPh™ PhJIoo ^ (271)^+-? ' 

so 

n 

where the constant Ck depends only on k (roughly, it equals the number of summands 
times k^). Hence, when taking n — )■ 00 we will not have to take care about Bm,n- 

Let us have a look at ([8]) and compute now the limit of the first term when — )■ 00. 
We observe that Iw,, ~^ 1 pointwise on M*'. Moreover, ly^ — )■ l[o,27r)P5 and ly^ — )■ for 
r] such that rjs 7^ r]t but 7r(s) = 7r(t) for some s ^ t. Thus we consider only r^'s such that 
fls = Vt whenever 7r(s) = 7r(t) and then the following simple observation 



yields for all these 77's 

(-) 

\mn J 



'2tt 
mn 
p 

JJdet 

i=i 



27r 



mn 



{z, - zt) + 2iT{r]s - r^t) + - 0:^(5) 



i,t=l 



~ 



1 ^ 

in 



By estimate ([7]), PE^y) is bounded by k^^^/m^, so the integrand in ([8]) can be simply 

bounded. Thus by Lebesgue's dominated convergence theorem 



For any p < m the integral /[q 27r)p J^^^ equals m!/(m — p)!. Consequently, we 
finally obtain 

p 



p[Ji Jzi, . . . , 2;^) y Y] — 



m! 



ri-s>oo ■< — ' m' 
p,7r 



(m — p)\ 

^ ' 1=1 
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In view of ([5]) this completes the proof. □ 



2.2. Proof of Theorem ([b]). Fix a point z = {zi, . . . , z^) E M'^. We let m and n tend 
to infinity and want to prove that p^^ (z) tends to 1. Recall ([HD and notice that due to 
estimate ([9]) all the terms with p < k — 1 are bounded above by Ck/m, so we can write 

k 

- ) Psli^)p3' iyiz,x))dx. 



Using the formulas for the correlation functions and the permutational definition of the 
determinant, we can put the integrand in the following form 



27r 



1 ^ 




k 


■ det 






- s,t=l 


n 


s,t=l 



{271)^ 



27r 



[Xi - X^ii)) ■ —SniVi - yr{i)) 



1+ ^ sgncTSgnrJJ — s^(a 

a^id or r^id i=l 

where the second summation runs through permutations a and r of k indices. The point 
is that each term in this sum tends to zero with m and n going to infinity as we have 



27r 



SmX Xi 



— -7- for i such that i 7^ cr(i), and —Sniyi — yrii)) — — > if i 7^ r(z) 
(see ffTOj) and mind the fact that actually y depends on m and ra). Recall also that Iwr, ~^ 1 
and ly^ — 7- l[o,27r)'=- Moreover, ([6]) yields that the whole sum is bounded by {k\Y/ (27r)^. 
Therefore by Lebsegue's dominated convergence theorem we conclude that 



m.n-^co 



> / l[0,27r)'=(a;) 



(2vr)' 



-dx 



which finishes the proof. □ 



3. Concluding remarks 

At the very end we shall discuss the tensor product of more than two matrices. We only 
briefly sketch what can be easily inferred looking at the proof of the main result. 

Let El, E'^, E'^ be the point processes of eigenphases of independent I x I, m x m, and 
n X n random unitary matrices respectively. Proceeding along the same lines as in the 
proof of Theorem (jaj), we conclude that the point process ® E'^ ® — vr) locally 

behaves as the Poisson point process on M when / is fixed but m and n tend to infinity. 
Indeed, the asymptotics of the A;-th correlation function p^^^ (z) of that process is governed 
by the integrals 

/ \T~') Pa (a;)pg^ (?/)pP(w(x,?/,2;))da:d?/, 

i[o,2^]p+*nv, \lmnj 

which we then sum suitably. Expanding the determinantal correlation functions of E'^ and 
E'^ (see the proof of Theorem (jb])) we find that the limit of p'^^\z) equals J^jp^. ~ ^' 
where the last identity is due to the well-known combinatorial fact that Yl^p=i '^&{k,p)x{x — 
1) ■ . . . [x — p + 1) = x^ . The same line of reasoning applies also when in addition / — > oo. 
Then the asymptotics depends only on the integral 

/ (t"") Ph?(^)Ph^^ (?/)pS'(^(a;.Z/.^))da;dt/. 
J[o,2^]2fc \lmnj ' 

Again, we carry on as in the proof of Theorem (jbj. 

Let An], i = 1, , 2, . . ., be independent Ui x rii random unitary matrices. The other cases 
of tensor products ^^Li^nl, when for instance all but one of rij's are fixed, seem to be 
more delicate and we do not wish to go into detail here. Moreover, it looks challenging to 
consider the tensor products when the number of terms M tends to infinity and (^i)^i is 
fixed. The simplest case of rij = 2, i > 1 has been addressed in [TSKZZj . 
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